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ABSTRACT 


The radiative transfer of Lya, LyB, and Ha in a hydrogen gas 
containing dust is considered with application to the nuclear gas 
in Seyfert galaxies. The line radiation originates from the 
kinetic energy of suprathermal particles, lost initially to a free 
electron gas, which then delivers it to the radiation through 
excitation of neutral hydrogen atoms. By neglecting the direct 
escape of line radiation and by averaging over the gas , we suppress 
the transfer in space and consider transfer in frequency only. 

The dust degrades the line radiation via frequency-independent 
absorption, converting the energy to infrared luminosity. We 
solve for the source functions in the lines, using appropriate 
approximations, in order to determine under what conditions the 
narrow component of the Balmer line radiation from the gas can be 
self-absorbed and degraded without similar degradation of the broad 
component, which originates from the suprathermals themselves. The 
results of the transfer calculation are used to find self-consistent 
values for the temperature and ionization of the gas for various 
amounts of dust and various concentrations of suprathermal particles. 
Finally, we examine the luminosity in the Balmer continuum that 
would result under these conditions. 


I . Introduction 


If the broad component of permitted emission lines observed 
in quasars and Seyfert galaxy nuclei arises from suprathermal ions 
(Ptak and Stoner 1973, Stoner et al 1974), the ambient gas clouds 
in which these ions are slowing must be optically thick in the 
Lyman and Balmer series. This requirement follows from the necessity 
to trap and ultimately to degrade the emission from ambient, thermal 
atoms that would otherwise produce very strong central narrow cores 
in the observed line profiles. 

In a previous paper (Ptak and Stoner 1975) we discussed the 
physical conditions in an ambient gas penetrated by suprathermal 
protons, assuming all line emission from the ambient atoms is 
trapped. In order to further investigate the feasibility of the 
suprathermal protons idea, we deal here with the question of how 
to produce the necessary large optical thickness in the Balmer 
lines. Large infrared luminosities have been observed for several 
Seyfert nuclei and this suggests that the nuclear gas may contain 
dust. Therefore, in this paper, we consider cases where the bulk 
of the kinetic energy of the suprathermal particles is ultimately 
transformed to infrared luminosity by dust. 

The presence of dust is related to the problem of producing 
a sufficient optical thickness in Balmer lines since the n = 2 
population of the ambient hydrogen atoms depends on the trapped 
Lyman alpha intensity, and absorption by dust is a mechanism for 
d.cs'troyi.n^ Lynicin slph^* So v/s consider tho r 2 .di 3 .t 1 v 6 transfer of 
Lya (in frequency space) when a dust-like absorption mechanism is 
operating . 
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We present the results of this calculation as a function of 
the density and temperature of the electron gas and of a parameter 
which describes the amount of dust present relative to the amount 
of gas. 

A large optical depth is not sufficient, by itself, to insure 
that Ha photons are degraded before they leak out of the gas. In 
order to deal with the Ha and LyS radiation, we have converted the 
transfer equations for a three-level atom to an approximately 
equivalent two-level transfer equation, using the results of the 
Lya calculation to fix the n = 2 population. This allows us to 
examine the effect of the dust on the Ha source function. 

With the Lya source function detesniined, we can calculate the 
temperature and ionization of the gas in a self-consistent way, 
assuming that the cooling is dominated by infrared radiation from 
the dust. We have done the self-consistent calculation, considering 
the gas to be everywhere the same, and the results are presented 
for various amounts of dust. 

Since a small amount of continuum cooling also occurs in the 
cases we have considered, we discuss what the Balmer continuum 
emission should be like if the gas is thin to this radiation. 


Ppo.T)UCIBn.ITY OP THE 
Sg^i. pagb is poor 
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II. The Lya Source Function 

We wish to calculate the source function which describes the 
trapped Lya radiation field in a dense, optically thick gas that 
is heated by a flux of suprathermal particles. The density and 
the linear dimensions of the gas cloud are assumed large enough 
that the overwhelming fraction of the gas is at least several 
thermal ization lengths from the surface. In this limit the source 
function and the mean intensity of line radiation are equal nearly 
everywhere, so we can consider the source function to depend only 
on frequency. 

Under these conditions the radiative transfer problem for a 
two-level atom reduces to the condition for statistical equilibrium 
of the photons at each point in frequency space near the line 
center. The photons are assumed injected with a Voigt profile, 
then to migrate in frequency due to collisions with atoms until 
they are destroyed in one of three ways: collisional de-excitation, 

an absorption by an atom leading to the production of a 2-photon 
continuum, or absorption by dust. The profile for absorption by 
atoms is the same as the emission profile, while the profile for 
absorption by dust is assumed frequency independent. The redistri- 
bution in scatterings by atoms is assumed coherent in the atom's 
frame, the situation originally called Case II by Hummer (1962) . 

There are three parameters which determine the nature of the 
radiation field. The collisional excitation and de-excitation 
process is described by which is equal to the rate of collisional 

de-excitation normalized to the rate of Lya absorption by atoms. 
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The rate of two photon emission from the 2s state/ normalised in 
the same way, is the parameter effect of the dust is 1 

characterized by e^/ which is the ratio of the coefficient for f 

absorption by dust to the total coefficient for absorption by atoms 
at the line center, multiplied by <p Co) . | 

Equating the rate of absorption emd emission processes at | 

each frequency and at each point in space, the transfer equation | 

reduces to the form: 


S(x) = e3^4>(x) + 


[ 1 - - ej] r 


S(x')R(X/X')dx' (1) 


where 


V - V. 




AYjj being the Doppler width, S Cx) is the source function in units 
of the Planck level at the electron temperature, and (j)(x) is the 
is the normalized Viogt profile. For the angle-average redistribution 
function R(x,x'), we choose (see Adams et al (1971)). The 

A Q 

electron temperature and density determine at T = 10 K, 

= 2.1 X 10 For the densities, we are considering, 

collisional mixing of the 2s and 2p states proceeds very rapidly, 


and ^2 simply the ratio of the two-photon transition probability 

-9 

to that for 2p -► Is, €2 * 5 x 10 

We have solved equation (1) numerically, using the method of 
cubic splines suggested by Adams et al (1971) , choosing various 
values for N^, T, and e^. The source function is displayed in 
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figure 1 for four values of with the electron density and tem- 

9 4o 

perature taken to be 5 x 10 cm ” and 10 K. We find from our 
results that it is possible to characterize the source function 
with two simple, semi-empirical formulas for §^ 3 ^ 2 ' source function 
averaged over the line profile, and x', the "effective half width". 
This is the halfwidth of a hypothetical trapezoidal source func- 
tion with the same total area and altitude as the computed one. 

To a good approximation, x* is the halfwidth at half maximum of 
the source function. The semi -empirical fomulas are; 

* .15 £ 3 ^ , and (2) 

"1 

£l + £3 + 2x'£3 



( 3 ) 
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III. The Three Level Atom 

If we wish to Investigate the trapped Ha field, we must 
include the n = 3 states for the hydrogen atoms of the gas, so 
the LyB radiation field is immediately involved. Absorption of 
a LyB photon by a gas atom may be followed by an Ha emission, and 
likewise an Ha absorption may be followed by emission of LyB. At 
equilibrium the rates of these two processes will be the same; if 
we assume that this balance applies for each set of corresponding 
frequencies, we can simplify the three-level transfer problem. 

As in ihe case of Lya, it is convenient to define 
S 23 (x), the source functions for Ha and LyB, respectively, each in 
units of the Planck function at that frequency. The doppler fre- 
quency AVjj appropriate for each line is proportional to the line 
center frequency, so that an atom in the n = 3 state can emit 
either a LyB or an Ha photon in a given direction with the same 
dimensionless difference frequency x. A similar statement holds 
for absorption processes. Using the approximation that the n = 2 
population is independent of the speed of the atoms and is deter- 
mined by the Lyman alpha intensity alone, we find from the balance 
of Ha and LyB processes at each x that 

SisCx) = S^2 S23 (x) . (4) 

Balancing the absorption and emission of Ha and using equation 
(4) , we obtain the transfer equation which determines the frequency 
dependence of the Ha source function: 

RFP"’' 


i 
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S23 




e,*(x) + (1 


-.i.n 


833 Cx)R(x,x')dx' (5) 


The notation Is similar to that used In equation (1). Here R(x,x') 
Is Rjj ^ with an effective value of the parameter a, which Is ob- 
tained uslnq the lifetimes of both the n » 2 and n » 3 levels, 
appropriately weighted according to branching ratios, 


= 1.12 X 10“^ , 


The three pareuneters In equation (5) can be related to the 

parameters In the Lyman alpha transfer problem and to the resultant 

Lyman alpha source function by comoaring the rates of corresponding 

processes. We will write the relationships for these parameters 

4o 

for an electron temperature of 10 K, so the slight temperature 
dependence is not included in these formulas. Comparing the rate 
at which the electron gas will collisionally excite the n - 3 level 
with the rate of collisional excitation of n 2 and adopting 
2.0 X 10 for the ratio find that to two significant 

figures. 


ej^(Ha) « 150 {1 + 2.6 x 10"^/Sj^2^ e^^aya) 


( 6 ) 


Comparing similar collisional de-excitation rates. 


cJ^(Ha) * 150 ej^(Lo) (7) 

The dependence on ^12 in equation 6 derives from the direct exci- 
tation of n s» 3 from n » 1, which becomes relatively more important 


than the excitation from n ■ 2 when the Lya source function Is 
small. 
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If the dust particles are relatively large so that they 
present the same cross sectional areas for absorption of Lyman 
and Balmer radiation^ then a comparison of the absorption rates 
gives 

e^CHa) ■ 4300 S^2 * 63 (L«) , (8) 

where the dependence on the Lya source function results from 
choosing to measure the source functions In units of the Planck 
function at the line centers. Equation (8) should be regarded 
as an upper limit to the appropriate for equation 5/ since If 
the dust Is very small/ Its cross section for Ha relative to Lya 
could be reduced by a factor as small as the ratio of wavelengths. 
We have performed calculations using both the big dust and small 
dust limits. 

The actual calculation proceeds by choosing values for the 
electron density, the temperature and the dust parameter for Lya. 
Equations 2 and 3 then yield the properties of the Lya source 
function, which, together with equations 6, 7 and 8, give the 
parameters In equation 5. Equation 5 Is then solved by the method 
previously used to solve equation 1. 

The importance of the results of this calculation Is that 
they can be used to test the feasibility of degrading the "narrow 
core" radiation without doing the same to the broad line radi^4tion 
that is produced directly from the suprathermal atoms. We wish to 
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sett If the gas can be thin to dust yet thick to Ha at all frequencies 
where ' j Ha source function is large. So we have computed the 
ratio of the absorption coefficients of the at^s and the dust, 
t(x')/Tjj - ♦(xM/Cj at the "effective half width" frequency for 
S 23 . The ratio tends to be slightly greater than one when the 
effective half width is smaller than about 3.5 Doppler units. When 
the effective half width is greater than 3.5, a typical photon 
must undergo several nearly coherent scatterings before reaching 
this frequency and this results in higher values of the ratio 
t{x*)A|j» about 40 at x' ■ 10. 

Larger values of x' are most easily produced at higher 3 lec- 

tron densities, as the typical values in Table 1 show. Since 

4o 

x* » 10 corresponds to 100 km/sec at 10 K and since the HWKM of 

the "broad component" from supra thermal atoms is about 3500 km/sec, 

there is a large range of physical conditions where the dust can 

very effectively remove the "narrs^w core" without appreciably 

attenuating the broad component. 

We summarize the results of the calculations in table 1, 

where we display values of: the effective half-width of the Ha 

source function, the optical depth in Ha at this frequency, and 

the optical depths in the Balmer and Paschen continuums. These 

are all expressed relative to the optical depth in dust, and the 

4o 

gas temperature is taken to be 10 K. 

If the size of the region is just the distance needed for the 
suprathermal protons to stop in the gas, the neutral hydrogen 
column density can be simply related to the ionization of the gas 
and the initial energy of the most energetic supra thermals. From 
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20 

Ptak and Stoner (1973) / we have N^L « 10 
optical depth at the line center is 



Tq (L ya) 


« 10 ' 



and so the Lya 


O 

Taking as a typical value T^(Lya) « 5 x 10 , the values for 
can be converted into optical depths for the dust. We could 
then convert the optical depths in table 1 into absolute values. 
For this choice of T^(Lya)/ an of 10"^^ corresponds to a 

• Q 

of about 0.1 f ■ 10 corresponds to Tq « 1.0, and so it goes. 
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IV. Self -Consistent Physical Conditions 

Having obtained the depetidence of the Lya source function on 
the electron temperature and density and on the dust parameter 
we can now determine self-consistent values for the temperature 
and ionization of the gas by solving the appropriate balance 
equations. In doing this^ we assume that the gas is heated by the 
suprathermal particles which are thermalized in it. and that the 
radiation field is determined by the suprathermals and by the gas; 
there are no external sources. 

In photo-ionizaticn models (MacAlpine 1972) . a large flux of 
ionizing photons is present in the line-emitting region of quasars 
and Seyfer>* galaxies. The existence of a. hard uv flux does not 
follow directly from observations, but it is assumed to exist from 
extrapolation of observed optical continue and in analogy with 
other, more familiar objects such as planetary nebulae. If there 
is such a photon flux in any of these objects, then we must con- 
sider what effect it would have on the calculation we are doing 
here. Let us consider a dense cloud of gas at some distance from a 
central continuum source. Incident on one side of the cloud is 
the ionizing photon flux and the flux of suprathermal particles. 

If the gas density is large enough, the particles will penetrate 
further than the photons, and it is in the region shielded from 
the photons that our calculation applies. 

Let us make a numerical estimate of the density required. 

The photo-ionized layer will have a thickness d determined by; 

P - d . 


(9) 
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2 

where F is the number of ionizing photons incident per cm per 
sec., Ng is ‘-he electron number density, and is the total re- 
combination coefficient. The suprathermal protons will penetrate 
a distance d* given by 

d'fe 10^° ^ , (10) 

6 

where is the initial kinetic energy of the protons (see Ptak 

and Stoner 1973) . Let us suppose that the cloud is one light year 

from the source, then, for an ionizing photon luminosity of lO^^sec”^, 

19 —2 -1 

typical of a bright quasar, F « 10 cm sec . If we take 

10 -3 

10 MeV, then d' exceeds d if > 10 cm . For a typical 

IQ —2 —1 

Seyfert galaxy Ituninosity, F « 10 cm sec , and we need 

9 —3 8 —3 

Ng > 10 cm . These densities are within the range of 10 cm to 

11 -3 

10 cm which we are considering in this model, and so the 
assumption of no external source is a reasonable one. Of course, 
a more complicated geometry could reduce the required density by 
a large amount. Also, It is possible that the optical continuum 
arises from a hot gas which is energized by the suprathermal 
particles themselves, and which is reddened by the dust. 

We find the temperature and ionization of the gas by simul- 
taneously satisfying two balance equations. The ionization balance 
is described by equation (1) of Ptak and Stoner (1975) , except 
now the Lya number density is calculated using the correct source 
function. The energy balance equation is obtained by equating the 
rate at which the gas is heated by the suprathermals with the rate 
at which it is cooled by the dust; 


UCUiii,,; ; -/o' llii- 
OKIGiNAL PAGi: . iOin; 
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NeE„ 

— - 0 c N It (11) 

T 

Here Ng is the number density of suprathermals/ t is the time for 
a suprathermal to stop in the gas, U is the Lya energy density, 

Nq is the dust number density, and 9^ is the dust absorption cross 
section. If we take the initial energy to be about 1 MeV, this 
equation becomes 

Ng f^(2 X 10“^^) = (.031) (2 • X* • - f^) 

( 12 ) 

X Sj exp(-Ej^2/^^i » 

where f^ is the ionized fraction of the gas, and 10.2 eV. 

Solving these balance equations together with the use of 

equations (4) and (5) gives us the ionization and temperature of 

the gas for given values of N-, N„, and e . Representative re- 

on 3 

suits are presented in table 2. 

For those cases which are displayed, a small amount of the 

cooling is accomplished by continuum processes (a few percent of 

the total) . Therefore, if the gas is i:hin to the Balmer continuum, 

an appreciable amount of energy comes out in this way, otherwise 

it goes into the Paschen continuum. We can estimate the relative 

strength of the Balmer continuum to that of H0 fairly simply. The 

emission rate of Balmer continuvim photons is about the rate of 

2 

recombination to the n = 2 state; a~. In the context of our 

e 2 

model , the HB emission s due to the suprathermal protons and the 
NR 

rate is =-, where Rg is the number of H6 photons produced by 

T 
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each suprathermal which stops In the gas. For the small ionization 
fractions indicated in table 2, Rg * 10. If we divide these two 
rates, we find that the number of Balmer continuum photons pro- 
duced per R6 should be 


L(BAC) _ ^e^a"^ 

Lmr ‘ rrsg 


8 


"s i 


( 10 '^) 


(13) 


for Eq * 1 MeV. 


And since 


6 


this ratio must be on the order of 100 or more when f^ = 0.1. 

A large value for this ratio is apparently observed is some 
low redshift quasars by Baldwin (1975) , but the profiles observed 
for these objects are not fit by profiles we have calculated so 
far. We know of no published values for this ratio for those 
Seyfert galaxies and quasars for which we have obtained a com- 
pelling fit (for these fits, see Ptak and Stoner 1973 and Stoner 
et al 1974) . 
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V. Discussion of the Results 

We have investigated the conditions under which a dense gas 
heated by suprathermal particles and cooled by dust will be op- . 
tically thick in L}^n and Balmer lines. At the same time, we 
have determined values of the temperature emd ionization of the 
gas that are consistent with such a picture. These two aspects 
eure intimately Interconnected; a large value for the dust para- 
meter will produce a fairly neutral gas for a wide range of number 
densities for the gas and the suprathermal particles. On the 
other hand, a large lowers the optical depth in Ha relative to 
that of the dust. 

In order for the cooling of the gas to be accomplished by 

infrared radiation from dust, the parameter cannot be much less 

than 10*"^^. The need to trap Ha photons until they can be degraded 

—8 

requires that be about 10 or smaller. 

The calculations reported here appear to establish the feasi- 
bility of producing broad emission lines like the ones observed in 
the nuclei of some Seyfert galaxies and in some quasars via the 
interaction of suprathermal atoms with a dense gas if an appropriate 
amount of dust is present to eliminate the "narrow core" photons. 
While these calculations have taken detailed account of the migra- 
tion of the line radiation in frequency due to the interaction with 
a thermal, partially ionized hydrogen gas, some other important 
features of a more realistic radiative transfer calculation are 
missing. 

Since the broad line emitting regions of these objects are 
unresolved and since quite different geometries can give the same 
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profiles, we have no strong clues to guide us in extending the 
calculation to include transfer in space. A full radiative transfer 
calculation, including transfer in space, seems uninteresting in 
the absence of stronger evidence for a particular geometry. 

Further, the transfer in frequency itself may be somewhat 
different from what we have assumed. The neglect of leakage of 
narrow core photons from deep inside the gas cloud is justified 
by the results of the calculation, but a large fraction of those 
produced near to the gas surface must escape easily. More important 
effects may come from deviations from a Maxwell -Boltzmann velocity 
distribution due to the presence of suprathermal atoms. For 
example, the enhancement in the "tail" of the distribution repre- 
sented by suprathermals could change the rate at which photons 
migrate from the line center to the wings of the line. Also, the 
atoms in the high velocity tail could directly produce a contri- 
bution to the broad component via the process we called "optical 
reverberation" in a previous paper (Ptak and Stoner 1975) . The 
effect of the latter on the source function would be similar to 
the effect of dust we have included in the present calculation. 

The full range of radiative transfer effects due directly or 
indirectly to the presence of suprathermal particles are probably 
capable of explaining the large Balmer decrement observed in most 
broad-line emission objects. Such effects might also produce line 
profiles of quite different character from the asymmetric, flat- 
topped variety that our original simple models gave and thus be 
consistent with some of the more recent profile observations, so 
long as the widths of the observed profiles are not too large. 
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On the other hand/ it Is difficult at present to see how to 
produce hydrogen emission lines as broad as 15/000 km/sec or greater 
by the actions of suprathermal particles without invoking some kind 
of mass motion of the ambient gas (rotation/ high velocity expan- 
sion/ etc.)* This seems to be an undesirable increase in complexity 
in the model. However/ if the detailed profiles can be fit by 
assuming a single motion for all of the ambient gas or double 
sources with relative motion/ then the number of parameters would 
only be increased by one. 

We are in the process of including the effect on the radiative 
transfer problem of the suprathermal atoms in a self-consistent 
way. This can be important both for the energy balance of the gas 
and for the radiative transfer problem. It is possible that for 
those gas densities at which the suprathermal atoms act as the 
major cooling agent (the optical reverberation process) / the re- 
quired flux of suprathermal particles and the line profiles produced 
will both be substantially affected. The results of this self- 
consistent calculation will be presented in a subsequent paper. 
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FIGURE CAPTION 


The frequency dependence of the Lyman a source function for 

•7 

four different values of the parameter using = 10 and 
£2 * 5 X 10 . The source function Is in units of the Planck 

level so S " 1 corresponds to thermal equilibrium between the 
radiation and the free electron gas. 



FREQUENCY IN DOPPLER UNITS 




TABLE 1 


Optical Depths Relative To Dust 
















TABLE 2 


Temperature and Ionization of the Ambient Gas 
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®3 

T 

‘i 

lo’ 

10 

10“® 

12,650 

.31 

10® 

10 

10“® or 10“^° 

10,040 

.09 

10® 

10^ 

lO"® 

10,060 

.09 

10® 

lo2 

10-®® 

10,010 

.13 

io“ 

10® 

10~® 

10,310 

.13 

10® 

10 

10“® or 10“^° 

8,740 

.03 

10^® 

lo2 

10"® or lO"^® 

8,850 

.03 

io“ 

103 

io“® 

9,250 

.03 

10^1 

10^ 

10-1® 

9,940 

.07 

10® 

1 

10"® or 10"1® 

7,930 

.008 

io“ 

10 

10"® or 10"1® 

7,950 

.008 

10^1 

102 

10"® 

8,100 

.008 

io“ 

10^ 

10-1® 

8,390 

.01 

lol® 

10^ 

lO"® 

8,830 

.00? 

10^2 

10^ 

10-1® 

9,520 

.02 








